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ABSTRACT 

We study the /5, N critical behaviour of non compact QED with iV 
species of light fermions, using a method we have proposed for unquenched 
simulations. We find that there exist two phase transition lines: one, 
second order, and the other, first order, that approaches asymptotically 
the (3 = axis. These two lines have different physical origin, the second 
one being entirely due to fermions effects. We discuss the effect of the 
approximation used, in terms of an expansion of the effective action in 
powers of N, and conclude that the general features should not be affected 
by this approximation. 
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Strongly coupled lattice QED has been, in the last years, the object of 
much attention, both analytical and numerical. The main subject of these 
investigations has been the question if a strongly coupled abelian model 
can be constructed with non trivial interactions in the continuum limit. 

During some time, the results about the nature of the fixed point of 
this model were controversial in the sense that different groups reported 
different values of critical exponents respectively consistent or inconsistent 
with a mean field description of this model. More recently, an extensive 
investigation of the quenched model in large lattices and at small fermion 
masses [1,2] gave evidence that a precise determination of the critical 
coupling (3 C is needed in order to extract critical indices with good accuracy, 
since small variations in the value of f3 c may induce strong changes in the 
values of the critical exponents towards their mean field values [3]. 

The results of [1,2] reporting non mean field values for the critical 
exponents in the quenched model, have been improved for the unquenched 
case in [4,5] and also by us in [3,6] in a completely independent calculation 
which allowed the numerical determination of (3 C from the results of the 
plaquette energy without the need of any kind of extrapolation at zero 
fermion mass. Additionally, an analysis, based on Renormalization Group 
approach, performed in [6] on the results for the fermionic effective action, 
pointed again to a non gaussian nature for the strongly coupled fixed point 
of this model. These results are in contrast with those in [7] . 

The agreement between the results of refs. [4,5] and those of refs. [3,6] 
for the two and four flavour cases was certainly encouraging and suggested 
the existence of a non trivial continuum limit for strongly coupled QED, at 
least if the number of dynamical flavours is less or equal than four. 

Another intriguing result reported in [5] was that critical indices for 
the two and four flavours models, obtained from numerical simulations and 
the use of the equation of state, were found compatible between them and 
also with the critical exponents of the monopole percolation transition, thus 
suggesting that different flavour models are in the same universality class. 

This result can be also understood in our approach, if the critical 
value of the plaquette energy is independent on the number of flavours, as 
suggested by the results of [3]. In such the value (for instance) of 

the 8 exponent which measures the response of the system to an external 
symmetry breaking field will be given by the dominant contribution to the 
expansion of the chiral condensate in powers of the number of flavours N 
at fixed pure gauge energy E c [3] , and it should be independent of N. 

The aim of this letter is to analyse in detail the dependence on the 
number of flavours A of the physical results and in particular to investigate 
the phase diagramm of this model in the A, (3 plane with special attention to 
the analysis of the physical origin of the different phase transitions observed. 

The approach we use to simulate noncompact QED with dynamical 
fermions is that of refs. [3,6], originally tested in the compact model [8], and 
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based on the introduction of an effective fermionic action S E ^(E, N,m) 
which depends on the pure gauge energy E, fermion mass m and number 
of flavours N, and which is related to the gauge fields by the relation 

[3] 



e -S* ff (E,m) = < ( det A )T >E = 

t.(x)-WE) 

(1) 



JidA^xWetA^A^x^S^Z^F^-GVE) 



J[dAp{x)]6{i^ <v F*,(x)-6VE) 

where A(m, A^(x)) is the fermionic matrix (we use staggered fermions), 
E is the normalized pure gauge energy and the denominator in (1) is the 
density of states which can be analitically computed 

N(E) = C G E*<y-V (2) 

Cq in (2) is some irrelevant divergent constant and V is the lattice volume. 

After the definition of the effective fermionic action, the partition 
function of this model can be written as a one-dimensional integral 

Z = J dEN(E)e- 6f3VE - s "ff ( - E ' N ^ (3) 
from which we can define an effective full action per unit volume as 

S eff (E, 0, N, m) = - ^ In E + 6(3E + Sf fJ (E, N, m) (A) 

S E jf(E, N, m) in (4) is the effective fermionic action (1) normalized to the 
lattice volume. 

The thermodynamics of this system can be studied now by means of the 
saddle point technique. The mean plaquette energy < E p >= Eo(m, (3, N) 
will be given by the solution of the saddle point equation [3] 



AE f 6dE 
satisfying the minimum condition 



TB-P-e»*S?ff&,N,m) = (5) 



^ + l^ ff (E,N,m)>0 (6) 
By differentiating equation (5) respect to (3 we get for the specific heat 



C,= l p <E p >=-{ 



AE*(m,(3,N) 



+ l^2 § eff(E,N,m) 



E (m,(3),N 



} _1 

(7) 
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The effective fermionic action S^ff was computed in [3,6] as a power 
expansion on the flavour number N. 

The observation of two different regimes in the behaviour of as 
a function of E, linear in the small energy region and nonlinear at large 
energies, and a carefull analysis of the numerical results allowed us the 
determination of the critical value of the energy (which turned out to be 
independent on the number of flavours N) and of the critical coupling (3 C . 
Furthermore it was established in [3] that the results for Sfff reported in 
Fig.l can be very well fitted by two polynomials with a gap in the second 
energy derivative of S^ff at E = E c which manifests itself in the specific 
heat Op as a second order phase transition. 

The origin of this non analyticity in the effective fermionic action is not 
clear at present. However, our numerical results suggest that this singular 
behaviour could have the same origin as the monopole percolation transition 
of the quenched model extensively analysed in [9]. In fact if we take our 
value of the critical energy E c = 1.016(10) (independent on the number of 
flavours) and consider the zero flavour limit, we get f3 c = 0.246(2) in very 
good agreement with the critical f3 c obtained in [9] from the results for the 
monopole susceptibility in the quenched model. 

Going back again to expression (7) it should be noticed that a non 
analyticity of the effective fermionic action is not the only way to get a 
discontinuity in the specific heat. A discontinuity in Cp can in fact be 
produced by a zero in the denominator of (7) through exact cancellation of 
the two terms in this expression [3]. To this end, a negative value of the 
second energy derivative of the effective fermionic action is necessary. But 
this is just what happens in the large energy region as it can be deduced 
from the results for the effective fermionic action reported in Fig.l. 

The above arguments can be made quantitative and at the same time 
the critical number of flavours N c can be simply estimated. We start from 
the cumulant expansion of the effective fermionic action [3] 

-S* ff (E,N,m) = j < Indet A(m, A^x)) > E 
N 2 

H {< (lndet A) 2 > E - < lndet A >|} + ... (8) 

32 

where < O >e means the mean value of the operator 0(A tl (x)) computed 
with the probability distribution [dA^(x)}5(± Y. x ,^ <v F l»{ x ) ~ 6VE). The 
numerical results for the succesive terms in the expansion (8) in a 8 4 lattice 
and massless fermions have been reported in [3], and are here extended to 
larger energies. As a first approximation, we will consider only the first 
contribution to (8), which is reported in Fig. 1 for A = 4. 

Fitting the results for the mean logarithm of the fermionic determinant 
at m = by two polynomials: first order for E < E c = 1.016 and fifth for 
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E > E c we get for the second energy derivative of the effective fermionic 
action 



= 



(E < E c ) 



(10) 



d 2 sF _ N 



[-1.669 + 3.877E - 2A9AE 2 + 0.494E 3 ] (E c < E < 2.5) 



dE 2 eff 4 



A simple analysis of these results tell us that in order to compensate 
the pure gauge contribution to the denominator of the specific heat in (7) 
we need N = 13.1. If AT is large but less than 13.1, the height of the peak 
increases with iV and just at this critical value, the specific heat diverges. 
Now, if we increase again the value of N, there will be an energy interval 
where the denominator of the specific heat in (7) will be negative and 
therefore no solutions of the saddle point equations (5), (6) will exist in 
this energy interval. This means that these energies will not be accessible 
to the system and hence, for N > N c a first order transition will appear. 

The phase diagram of massless noncompact QED in the AT, (3 plane 
which emerges from our results is plotted in Fig. 2. The continuous (broken) 
lines represent first (second) order phase transitions respectively. The 
end point of the first order phase transition line is a second order phase 
transition point with a divergent specific heat. The first order line ends 
at some finite (3 since J^S^f = for E < E c . On the other hand, the 
second order line merges, for large N, into the first order one since E c falls 
into the energy interval not accessible to the system, which widens as N 
increases. 

Of course the quantitative results of this simple analysis could change 
if we take into account higher order contributions to the effective fermionic 
action (8), but the physics behind these phase transitions can be understood 
at this simple level. What is independent on approximations is the fact that 
there is no first order phase transition at (3 = for any finite value of N. 

This result can be rigorously proved since the effective fermionic action 
Sgfj(E, m, N) is bounded for any finite N. Therefore, when (3 goes to zero, 
the effective action (4) has only one minimum at E = oo, whereas for any 
finite value of E the effective action is finite. Thus, at (3 = there is no 
first order transition for any finite N . 

A stronger result can be proved, under the very natural assumption 
(corroborated by the experimental data) that the effective fermionic action 
is monotonically increasing with E, namely that atA^^oo < E >= for 
all (3. This implies that, in this limit only one phase (Coulomb) does exist. 

In fact, since the value of the fermionic determinant is bounded, it 
follows that 



e 7V<logdetA> < e -S^ ff (E,m) < gJVlogA 



■ ?n a x 



(11) 
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which implies that the fermionic effective action diverges linearly with N. 
As a consequence, for N — > oo 



'e// 



and 



w - In A?" + NK(E, m) 

E=7i 2 



E=0 



(12) 



5, 



e// 



E=E 



NK(Eq, m) 



The monotonicity of S^j-(E) then implies that 



N— >oo 



lim 5 e// 



-5, 



e// 



= OO 



V£ ^ 



(13) 



(14) 



implying that E = is the absolute minimum of S e ff for any (3 ^ 0. 
This, together with the previous discussion, implies that the first order line 
reaches asymptotically the (3 = axis. 

To compare with published results on this subject, Kondo, Kikukava 
and Mino [10] found, within the Schwinger-Dyson approach, a continuous 
phase transition line which approches the (3 = axis asymptotically. The 
general discussion of the above paragraph excludes such a behaviour for 
large N. 

On the other hand Dagotto, Kocic and Kogut [11], in the framework 
of a numerical simulation, found evidence for a second order transition line 
wich becomes first order at large N, crossing the (3 = axis (N c ~ 30 
at (3 = 0). The probable origin of the disagreement between this result 
and ours is that the critical (3 at N = 30 is so small that it is difficult 
to distinguish it from zero in a numerical simulation where metastability 
signal can be observed also at (3 = [11]. 

From a physical point of view, these two phase transition lines have a 
different origin. The second order line, as pointed before, could have the 
same origin as the monopole percolation transition of the pure gauge model 
[5,9] . Indeed our results for the critical energy and their independence 
on the flavour number N favour this interpretation. The first order line 
is on the other hand produced by pure fermionic effects. When the 
number of dynamical fermions increases, the fermionic contribution to the 
denominator of the specific heat becomes more and more important and 
has the correct sign to cancel the pure gauge contribution. 

To finish let us say that these results have been also confirmed by our 
numerical simulations of this model. In these numerical simulations we have 
taken into account the first contribution to the effective fermionic action 
(8). In our opinion, the contribution of higher powers in N is not likely to 
change the qualitative behaviour depicted above. 

In fact, since the effective fermionic action diverges linearly with N, 
is bounded as E — > oo and is monotonically increasing, then, barring very 
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peculiar behaviours, it is convex for large E (corresponding to (3 — > 0). In 
such a case, for N large enough the two terms in (7) can always compensate, 
so a first order transition line will be in general present, at large N and small 
(3. In conclusion, although the numerical structure of the phase diagram 
might vary, it is very likely that the qualitative features remain unchanged. 

This work has been partly supported through a CICYT (Spain) - INFN 
(Italy) collaboration. 
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FIGURE CAPTIONS 

1) First contribution to the effective fermionic action (Equation 8) in a 
8 4 lattice, m = 0.0 and N = 4. Errors are smaller than symbols. 

2) Phase diagram in the ((3, N) plane at m = 0.0. 
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